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pec ARATIONS ; De Largtive Part_of Module 

HSHATAN = Standard H Float ng Point Arc Tangent 
MTHSHATAN2 = Standard H F gate. count Arctangent With 2 Arguments 

HATAN RB = Special wit rout “ 

MTHSHATAND - Standard H Floating po nt Arc Tangent 
MTHSHATAND2 = Standard H Floating Point Arctangent With 2 Arguments 
MTHSHATAND_R8 = Special HATAND routine 
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IDENT /2-006/ e: MTHHATAN.MAR EDIT: RNH2006 


MARA AAABAAAAAAAAALASALALALELELES SASSER SESE REESE REE RRS ES ERASER EE EEC ERECT CEES 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
pORPORAT ihe NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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ofACeatys MATH LIBRARY 
ABSTRACT: 


MTHSHATAN is a function which returns the H floating point era tongent 
value (in radians) of its H floating point argument. MTHSHATAN2 is two 
araunent H floating arctangent. The call is standard coh enroreterence. 
MTRSHATAN_R8 is a Special routine which is the same as MTHSHATAN except 
that a faster non-standard JSB cali is used with the argument in RO-R3 
and no registers are saved. 


MTHSHATAND is a_ function which returns the H floating point arctangent 
value (in cecrees) of its H floating point argument. MTHSHATAND2 is two 
acguaent H floating arctangent. The call is standard Cat Leny-cererence. 
MTHSHATAND_R8 is a special routine which is the same as MTHSHATAND except 
that a faster non-standard JSB call is used with the argument in RO- 
and no registers are saved. 


VERSION: 1 
HISTORY: 
AUTHOR: 
John A. Wheeler, 03-Oct-1979: Version 1 
MODIFIED BY: 
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Edit history for Version 1 of MTHSHATAN 


1-003 - Adapted from MTHS$GATAN version 1-001. JAW 03-0ct-1979. 
1- - Addea degree entry points. RNH 15-MAR=1981 


Edit history for Version 1 of MTHSHATAN 


- Mt - Use G* “atigshg k for externsis. SBL a4~Aug- -1981 
- - ppenge d MTH$$AB ya to MTH$$AB_ATAN 29-Sep-81 

4 = Change MTHSHATAND2 entry to MTHSHATAN2D in "aeuae to conform 

with original specific ation, RNH 05-Oct-81 


$-00> Chan ed ISTB to TSTH in M THSHATAND. RNH 05-Oct-81 
-006 = Un-did edit 2-004 to contend with PL/1, 
- Modi tied bea erqunens logic to avoid a microcode bug in the 
FPA. RNH 18-Dec-8 
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2006 BECLARATIONS > Declarative Part of Mods “CoSEPc1oad S1i2¢:4o EMTARTS eRe IATMGATaGemaR:1 28° chy. 2-0C 
; BS ~SBTTL DECLARATIONS ; Declarative Part of Module | 
0; 
00 4 ; INCLUDE FILES: MTHJACKET.MAR, MTHATAN.MAR 
308 $8 ; EXTERNAL SYMBOLS: 
0 95 -DSABL GBL 
880 38 eEXTRN MTHSK_INVARGMAT 
00 9 -EXTRN MTHSSSIGNAL : Signal SEVERE error 
0000 98 -EXTRN MTHSS$AB_ATAN_V ; Gobal table used by all Arctangent 
8 8 109 $ routines. Part of MTHATAN.MAR | 
000 101 ° | 
44 19 ; EQUATED SYMBOLS: | 
000041FC 4 104 ° ACMASK = “M<IV, R2, R3, R4, RS, R6, R7, RB> ; ENTRY register mask, | 
000 105 :; int ovf enabled 
0000 106 
0000 107 ; 
0000 108 ; MACROS: none 
0000 109 ; 
Bo08 119 ; PSECT DECLARATIONS: | 
00000000 118 .PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT | 
0000 «11 3; program section for math routines 
0000 114 ; 
0000 115 ; OWN STORAGE: none 
0000 116; 
0000 117 ; EXTERNALS: 
0000 «(118 ; 
0000 119 -EXTRN MTHSSSIGNAL ; Signal a severe error 
0000 120 -EXTRN MTHSK_INVARGMAT ; Invalid argument to math Library 
0000 = 121 -DSABL GBL ; No other externals allowed 
0000 \¢¢ 3 
44 \$7 ; CONSTANTS: 
00000000 00000000 00000000 0000C001 ais} 163 H_M1.0: .H_FLOATING -1.0 
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-SBTTL MTHSHATAN = Standard H Floating Point Arc Tangent 


+4 
; FUNCTIONAL DESCRIPTION: 
: HATAN = H floating point arctangent function 


; HATAN is computed using the following steps: 
1. If X > 4 _then 


Sete Ge Ge Ge Ge Se 


a. Let 


1/X. 
b. Conpute HATAN(W) = 


ree 

C. set RATANOO 
2. If =< 

a. Obtain xHI "' table lo 


b. Compute Z = (X = 
€. ey ner ® 


d. oeeete HATAN(XHI) high table 4 pe + HATAN(XHI) will have 


f X <0 th 


CALLING SEQUENCE: 


- Compute HATAN(E) e 
of degree 


= WeP(Wee2), where P is a polynomial of 


ook 
XHI)/(1 
Z*P(Z 


“at 


= Pie = HATAN(W) 


So eXHID. 
#2), where P is a polynomial of 


er bits, AN_XHI_HI, and the low 


2 HATAN. XHI_HI + (HATAN_XHI_LO + HATAN(Z)). 


en 
X + X 


*Q(X**#2), where Q is a polynomial 


en 
Compute Y = we he Xi) using steps 1 to 3. 
- Set HATAN(X) = -Y 


Arctangent.wh.v = MTHSHATAN(x.rh.r) 
-or- 
CALL MTHSHATAN(Arctangent.wh.r, x.rh.r) 


Because an H-floatin 
returned as the firs 


to the second argument 


Standard. 
INPUT PARAMETERS: 

LONG = 4 

x = 2 * LON 
IMPLICIT INPUTS: 


G 


OUTPUT PARAMETERS: 
VALUE: H floating point arctangent (in radians) of the argument 


Arctangent 


none 


1 * LONG 


resul 
ot oe 


t cannot be expressed in 64 bits, it is 
ent, with the input parameter displaced 
accordance with the Procedure Cal ing 


; Define longword multiplier 
; X is the argument 


; Arctangent is the result 


Page 
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NOUS WN OO OONOUES Win 


re] 
FP ot at at at ot ot ot ot ot a a 


J1 
oint Arc Tangent Fun , a 16-SEP=1 
andard H Floating Point Ar 6-SEP-1 


IMPLICIT OUTPUTS: none 
SIDE EFFECTS: 
Signals: none 
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; NOTE: This procedure disables floating point underflow, enables 


; integer overflow, causes no 


ENTRY MTHSHATAN, ACMASK 


MTHSFLAG_JACKET 


MOVAB G*MTHSSJACKET_HND, (FP) 


MOVH  ax(AP), RO 
BSBW § MTHSHATAN RB 


MOVO RO, @Arctangent (AP) 


RET 


floating overflow or other arithmetic 
; traps, and preserves enables across the call. 


andard call-by-reference entry 
sable DV (and FU), enable IV 
ag that this is a jacket procedure in 


t handler address to jacket 
ndler 


case of an error in special JSB routine 


Call special HATAN routine 
Store result in first argument 
Return to caller 


eee eee 


LT 
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Because an H floating result cannot be expressed in 64 bits, it is 
returned as the first argument, with the input parameter displaced 
to the second argument, in accordance with the Procedure Calling 


Rar sa ~SBTTL MTHSHATAN2 = Standard H Floating Point Arctangent With 2 Arguments 
aAr : FUNCTIONAL DESCRIPTION: 

gar ; HATAN2 =H floating point arctangent function 

Ar : HATAN2(X,¥) is computed as following: 

9A7 : If Y = 0 or X/Y > 2"8114, HATAN2(X,Y) = PI n xX) 

9A7 ; If Y > 0 and X/Y =< 2ee1{4, HATANG (x. ¥) = ) 

4 : If ¥ < 0 and X/Y =< 2*#114, HATANO(X,Y) = gn X) + HATAN(X/Y) 
9A7 : 

4 3 CALLING SEQUENCE: 

Ar : Arctangent2.wh.v = MTHSHATAN2(x.rh.r, y.rh.r) 

09A7 : -or- 

O9A7 : 

O9A7 : CALL MTHSHATAN2(Arctangent2.wh.r, x.rh.r, y.rh.r) 


OOOO OOO 00009000900 0909 09 0009 SIN NII NINN NINO OOOO OAOOOMIII «> 
NAUSEA © OD NOAU EWN OS OD NAUES WN 0 OONOUES WN —OOOnNOu OO 
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Q9A7 Standard. 
09A7 
Bone INPUT PARAMETERS: 
00000008 09A7 x = ¢ * LONG :; X is the first argument 
0000000C Seay y = 35 * LONG ; Y is the second argument 
Baar ; SIDE EFFECTS: See description of MTHSHATAN 
O9A7 > OUTPUT PARAMETERS: 
00000004 O9A7 " Arctangent2 = 1 * LONG : Arctangent2 is the result 
09A7 i-- 
Sn 
41FC O9A7 -ENTRY MTHSHATAN2, ACMASK ; Standard call-by-reference entry 
O9A9 ; Disable DV (and FU), enable IV 
S74 MTHSFLAG_JACKET ; Flag that this is a jacket procedure in 
6D Q0000000'GF 9E O9A9 MOVAB G*“MTHSSJACKET_HND, (FP) ; 
98 ; set handler address to jacket 
- ; handle 
98 498 3; Case of an error in special JSB routine 
50 08 BC 70FD 098 499 MOVH ax(AP), RO ; RO/R3 = arg! 
54 OC BC 70FD 1s #4 89 MOVH @y(AP), R4 ; R4/R? = arg2 
BA : : Test if Y = 0 or X/Y¥ > 20114 | 
42 13 O9BA ’ BEQL. INF ; Branch to INF if Y = 0 | 
58 50 $000 HH AB 098 5 BICW #*x8000, RO, R8 ; R&B = exponent (X) 
7E = 54 000 8F AB O9C 6 BICW #*x8000, R4, -(SP) 3; Stack = exponent(yY) 


| 
MTHSHATAN ;_H Floating Point Arc Tangen bok 16-SEP-1984 AX/VMS Macro V04-00 Page 14 | 
2-006 MTHSHATAN2 Standard H Pigating Point 8 orSEp= 1 984 9 3 9 MTMRTL. SREINTHHATAN. MAR; 1 . $63) 
8 8 A 9c8 7 SUBW2 (SP)+, R8 ; RB = ex ooeng (x) = exponent (Y) 
0073 af 1] af 9C 8 CMPW RB #115 : Com mpare R with 115. > | 
2c 14 44 5 BGTR 3 If xX/¥ > eniie, branch to INF 
20 1} : Test if Y>Oory<0 | 
54 BS 09D 18 , TSTW RG ; Test the sign of Y 
20 14 0904 214 BGTR A2PLUS ; Branch to A PLUS 14 Y>0 
0 B5 4 15 TSTW ; Test the s alga 8 of X 
OE 18 4 21§ BGEQ AIPLUS ; Branch to AIPLUS "e X >= 0 
+ 43 $18 : Y < 0 and X < 0 and X/Y =< 2##114 
4A 10 OQ9DA 520 $ BSBB MTHSHATAN_R8D ; RO/R3 = HATAN(X/Y) 
50 FAAF CF 62FD O9DC 521 SUBH2 4H P ; RO/R3 = -PI + HATAN(X/Y) 
04 BC 50 7DFD BaeS 5 $ MOVO RO, @Arctangent2(AP) ; Store result in first argument 
04 BoE : ? RET ; Return to caller 
0968 525 : Y <0 and X > 0 and X/Y =< 208114 Wo 
b9E8 3 $ A1PLUS: H_M 
3c 636-10 ~409E8 8 528 BSBB_- ss MTH$HATAN_R8D 3 RO/R3 = HATAN(X/Y) HM 
50 FAA CF 60FD O9EA 529 ADDH2 HPI, ; RO/R3 = PI + HATAN(X/Y) He 
04 BC 50 7DFD O9FO 530 MOVO RO, @Arctangent2(AP) ; Store result in first argument HP 
04 ber? 23) RET ; Return HP 
O9F6 338 : Y > 0 and X/Y =< 2#114 ay 
bors 33 A2pLus INF 
2E 1 O9F6 536 BSBB MTHSHATAN_R8D ; RO/R3 = HATAN(X/Y) 
04 BC 50 7DFD O9F8 537 MOVO RO, @Arctangent2(AP) ; Store result in first argument 
pore 238 ; RET ; Return 
O9FE 540 : Y = 0 or X/¥ > 2e#114 
SED ine 
BS O9FE 278 " -TStw RO : Test the sign of Xx 
OA 14 QA0DO 544 BGTR 1$ : Branch if X > 0 
10 13 pads 545 BEQL 2$ ; Branch if X = 0 
04 BC FAA7 CF 7DFD A04 546 MOovO H_MPI1_OVER_2, anrctangené2(aP) 
AOB 547 :; Store result in first argument 
04 ae 248 RET ; Return 
04 BC FA8F CF 7DFD QAO 239 1$: MOVO H_PI_OVER_2, BAre tengents(ar 
Ne 31 ; Store result in first argument 
04 OAl 23 RET ; Return 
Ale 88% se 
Aid He : Here if both X = 0 and Y = 0. Signal INVALID ARG TO MATH LIBRARY 
Ald = 557 
: #15, #1, R ; RO/R3 = reserved operand, copied 
ettne - sg At $58 is 9 . . ; to CHFS_MCH_S SAVRO/RI so handlers 
, ne 3 fae R2 3 can change Tf they want to continue. 
C OA L 3 
7 ‘BF QA OQAIA 6 MOVZBL rt hdae + INVARGMAT, -(SP) : Code for invalid argument to math Library 
00000000'°GF i FB OATE 88 CALLS , G*ATHS$SIGNAL 3 Signal SEVERE error 


THSHATAN Floating Point Arc Tan ent F ie tons 1 4 4: AX/VMS cro v04-00 Page 1 
" ATHSHA and 2 Standard H Ptgst ng Point g-8 P=] 9B4 iE 34 48 MTHRTL.S RCIMTHHATAN. MAR; 1 . ( 


04 OA25 564 RET ; Return if a handler says SS$_CONTINUE 
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. § $6 -SBTTL MTHSHATAN_R8 = Special HATAN routine 
" é 08 ; Special HATAN = used by the standard routine, and directly. 
A 53 ; CALLING SEQUENCE: 
A26 71 ; save anything needed in RO:R8 
A26 if 3 MOVH : Input in RO/R3 
Ace 60573: JSB MTHSHATA 
" ? Ae : return with we a RO/R3 | 
A26 278 ; Note: This routine is written to avoid causing any sateper overflows, floating | 
" e uf : ever tlaut, or floating underflows or divide by 0 conditions, whether enabled or 
OA26 579 ; 
OA26 580; mebahahd >): yee) 
_ 6 381 3 RO/R3 = floating argument then result 
A26 256 3 R4/R7 = scratch 
QOA26 583; R8 - Pointer into HATAN_TABLE 
OA26 38584 ; RO:RS = POLYH 
QA26 585 ; stack = Y during POLYH 
OA26 586 
QA26 = 587 
OA26 = 5588 
OA26 589 MTHSHATAN_R&D: ; For local use only! 
50 54 66FD OA26 590 VH R4, RO 
OA2A 591 MTHSHATAN_R8:: sé octal HATAN routine 
50 73D OAZA 592 TSTH RO | i s 1 @ arqment | 
03 18 OA2D 59 BGEQ POS_ARG : 
0087 31 pA 4 29¢ BRwW NEG_ARG ; Branch to negative argument logic 
pags 596 ; Argument is positive 
0A3 597 ; 
OA32 598 ®oS_ARG: 
58 50 03 9C QA32 599 ROT #3, RO, R8 i R8 = Biased exponent bits and first 
58 FFF80000 BF CA OA36 600 BICL  #@XFFF80000, RB ; | three fraction pits 
58 OOOIFFEC BF C2 OA3D 601 SUBL #*X1FFEC, RB ; Argument is less than 3/32, 
4C 19 QA4&4 o0¢ BLSS SMA : branch to small argument logic 
58 37 D1 OA4 60 CMPL #*x37, RB ; Argument is greater than of equal to 
4A 15 es one BLEQ LARGE _ARG $ 11, branch to large argument logic 
0A4B 606 ; Logic for positive medium sized arguments. Get pointer into HATAN_TABLE. | 
57 Q0000000'GF oD AS 608 MOVAL G*MTHSSAB_ATAN_V, R7 ; R7 = Address of RTL vector ntab 
57 saat bases *GF § AS 609 ADDL TY te ATAN_ V, R7 ; R7 = Address of MTHS$$AB Nm table 
6748 AS 610 MOVB (R7)CR Re 3; RB = Index into HATAN_TABL 
58 FSAD CF48 7EFD a2? el) MOVAO HATAN Robie R8), R8 ; RB = pointer to XHI 
Ab4 gig 3; Compute Z 
A6S 614; 
88 14 Ab4S 615 MOVO (RB) + 3 R4 = XHI 
4 63FD 0A68 616 SUBH3 RG, RO, Re SP) : (SP) = X = XHI | 
54 64FD OA6D 61 MULH ate RO ; RO = X*XHI 
98 60FD OA71 o13 ADDH RO ; RO = 1 + X*XHI 
66FD wr3 air DIVH +40 (SP) ; (SP) = 2 = (X = MHID/(1 + X®XH1) 
A79 631 | : Evaluate Z*P(Z*#2) 
A79_ 622 ; 


B 1 
MTHSHATAN ; H Floating Point Arc Tangent 4 & ES 16-SEP-19 1:34: 
2-006 ATHSHATAN.RB = Special HATAN routine 6-SEP- 198 9 3 98 
50 6E SFD OA79 623 MULH3 (SP), RO : 
F828 CF OE: $f sere OME 6 4 POLYH ty AHATANLENI=1, nATANTABT 
50 BE 64FD Age 626 MULH2 (SP)+, RO : RO 
20 88 60FD Ag9 627 ADDH (R8)+, RO : RO 
0 68 60FD OABD 628 ADDH (RB), RO : R 
A91 629 ; 
05 0A91 6 0 RSB : Re 
Oa 6 
A926 : 
OOA9 31 OA92 633 SMALL: BRW SMALL_ARG : Du 
OA 634 : 
A9S—s« 635 3 
bags B39 EL iti t test 
; Large positive argument Logic. 
OA9S§ §« 638 gs ; . 
0A95 = 639 
0A95 640 LARGE_ARG: 
F565 CF 50 67FD OA95 641 DIVH3 RO, H_M1.0, = (SP) i 6S 
50 6E 6E 65FD ASC 64g MULH C5) “(SP); 
F808 CF OE 50 75FD OAAT 643 POLYH RO, PHATANLENT- -1, HATANTAB| 
50 BE 64FD OAA8 645 MULH2 (SP)+, RO ; RO 
50 FAIF CF 60FD OAAC 646 ADDH2 H_PI_OVER_2_L0, RO : 
50 FAQS CF 60FD OAB2 647 ADDH2 4H ~PI“OVER~2-HI, RO : RO = 
05 OAB 648 RSB 3; Re 
OAB9 649 
OAS G81 f Logic f 
; Logic for negative arguments 
OQAB9 o3¢ ; 
OAB9 = 65 
OAB9 654 NEG_ARG: 

58 50 03 9C OAB9 655 ROTL #3, RO, RB : RB 
FFF80000 BF CA OABD 656 BICL AXFFFAOOOO, RB ; 
OOOSFFEC BF C2 OACS 657 SUBL #*XSFFEC, RB : Ar 

CS 619 OACB 658 BLSS SHALL : 
58 37 D1 OACD 659 CMPL =so#* RB : Ar 
48 15 OADO 660 BLEQ ON “TaRGES _ARG ; 
QAD2 661; 
m4 666 3 Logic for negative medium sized arguments. 
Q0000000'GF DE OAD2 664 MOVAL G*MTHSSAB_ATAN_V, R7 ~~; R7 
00000000'GF ¢O OAD9 665 ADDL GORTHSSAB ATAN-V, R7 —s;:-R7 
8 6748 90 OAEO 66 MOVE  (R7)CR RB : R 
58 F526 CF48 7EFD OAES 66 MOVAO NATANS ROASLE R8), RB : R 
OAEB 668 ; 
AEB 669 ; Compute Z 
AEB 670 ; 
54 88 70FD OAEB 671 MOVO. = (RB) + > RS 
7E 50 34 61FD OAEF 67¢ ADDH3 = R4, RO, Re SP) : (S$ 
50 «54 64FD (OARS) = 67 MULH2 RR : RO 
50 08 50 63FD OAF8 674 SUBH3 =oRO.. #1 : RO 
6 50 66FD Ar 679 DIVH2 RO. (sp)? > (S 
44 677 : Evaluate Z7*P(Z**2) 
50 6E€ 6€ 65FD 0801 ers MULH3 (SP), (SP), RO : RO 


—— see 
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= J**2 


P(Z*#2) 


HATAN(Z) = Z*P(Z**#2) 


turn 


mmy Label used to avoid adding 
an extra insrtuction in the 
medium argument logic 


P) = <W = -1/x 
= Weed 
= P(Wee2) 
= HATAN(W) = -WeP(Wee2) 
HATAN(X) = PI/2 = HATAN(W) 
turn 


= Biased exponent bits and first 
three fraction bits 
gument is less than 3/32, 
branch to small argument logic 
ouent is greater than of equal to 
branch to large argument logic 


Get index into HATAN_TABLE. 


Address of RTL vector entr 
Address of MTH$$AB_ATAN table 
Index into HATAN_TABLE. 
pointer to XHI 


HI = X = (=XHI) 


x 
*XHI = 1 + X*(=XHI) 
(X=(=KHI))/(01#X*(-xXHI)) 


~~ 
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POLYH RO, #HATANLENI-1, catia | 


MULH (SP)+, R 
(R8)+, R 
(RB), RO 


; Logic for large negative arguments 


RO, H_M1.0, =(SP) 
MULH3 (SP),~(SP). RO 


POLYH RO, #HATANLEN1=1, HATANT 


MULH2 (SP)+, RO 
SUBH2 H_PI_OVER_2_LO, RO 
SUBH2 H~PI-OVER@2-HI, RO 


> Small argument logic. 


SMALL_ARG: 
MOVO 


MOVO (SP)+, RO 


" 
bad 
=z 
eG 
, 
= 
~m 


Bete Se Be Se Fee te 


uM 


0 = HATAN(X) 
Return 


) 
Q = ix! 
Compare 2°-56 to :X: 
Needs polynomial evaluation 
Return with answer equal 
to argument 


RO, RO Pi RO = 
POLYH RO, #HATANLEN2-1, aa. 


MULH2 (SP), RO 
ADDH2 = (SP)+, RO 


n“~— 
oc 


T 

AN_XHI_LO + HATAN( 
HATANTX) = HATAN_XHI 
(HATAN_XHI_LO + RATAN 
; Return 


= HATAN(W) = PI/2 


2) 
$ HATAN(X) = KX + X#Q(X*##2) 
; Return 
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p67 d 1 -SBTTL WTHSHATAND = Standard H Floating Point Arc Tangent 
bee? 7? : 
p80! 724 34+ 
+U 4 ; ; FUNCTIONAL DESCRIPTION: 
dB? f 4 ; HATAND = H floating point arctangent function 
th 4 ; ; HATAND is computed using the following steps: 
0B67 731; 1. If X > 11 then | 
0B67 7 3 ; a. Let W = 1/X, 
0B67 733; b. Compute HATAND(W) = WeP(We*2), where P is a polynomial of 
0B67 734; egree 13. 
0867 735; Se ef HATAND(X) = pi/2 = HATAND(W) 
QOB67 736; 2. If 3/32 =< X =< 11 then 
0B67 4 =—737 ; a. Obtain XHI by table look-up. 
0B67 738; b. Compute Z = (X = XHI)/(1 + X*XHI). 
0B67 739; c. Compute HATAND(Z) = Z*P(Z**2), where P is a polynomial of 
0867 740 ; gegree 13. 
0B67 741 ; d. Obtain HATAND(XHI) by table look-up. HATAND(XHI) will have 
0B67 ee : two parts - the hioh order bits, HATAND_XHI_HI, and the low 
0B67 743 ; order bits, HATAND_XHI_LO. 
0B67 744 e. Compute HATAND(X) = HATAND_XHI_HI + (HATAND_XHI_LO + HATAND(Z)). 
0B67 8745 ; 3. If 0 =< X < 3/32 then 
0B67 746 ; a. Compute HATAND(X) = X + X#Q(X**2), where Q is a polynomial 
0B67 747 ; of degree 
0867 748: «= 4. IF X < 0 then 
0B67 38749 ; a. Compute Y = HATAND(!X!) using steps 1 to 3. | 
0867 750 ; b. Set HATAND(X) = -Y. 
0B67 751; 
QB67 ee ; 
0B67 753 ; CALLING SEQUENCE: 
0B67 754 ; 
0B67 755 ; Arctangent.wh.v = MTHSHATAND(x.rh.r) 
0B67 756; 
0B67 757; -or- 
0B67 758 ; 
444 $23 : CALL MTHSHATAND(Arctangent.wh.r, x.rh.r) | 
0B67 761 ; Because an H-floating result canmot be expressed in 64 bits, it is 
0867 bog : returned as the first argument, with the input parameter displaced 
0B67 765 ; to the second argument, in accordance with the Procedure Calling 
§867 764 ; Standard. 
67 765; 
0867 166 ; INPUT PARAMETERS: 
0B67 767 ; : 
00000004 0B67 748 LONG = 4 ; Define longword multiplier 
00000008 a4 168 x = 2 * LONG ; X is the argument | 
O86? 77 = IMPLICIT INPUTS: none | 
0B67 re ; 
au ae ; OUTPUT PARAMETERS: 
B67 a8 3 YALUE: H floating point arctangent (in degrees) of the argument 
00000004 67 in6 Arctangent = 1 * LONG : Arctangent is the result 
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B67 78 H 
$4 ag ; IMPLICIT OUTPUTS: none 
B67 781 ; SIDE EFFECTS: 
$867 ? : Signals: none 
0867 78s ; NOTE: This procedure disables floating point underflow, enables 
0B67 786 ; integer overflow, causes no floating overflow or other arithmetic 
0B67 787 ; traps, and preserves enables across the call. 
0B67 788; 
B6/ 44 jeer 
B6 790 
0B67 791 
41FC 0867 oe -ENTRY MTHSHATAND, ACMASK ; STANDdard call-by-reference entry 
0869 79 ; Disable DV (and FU), enable IV 
p86? 794 MTHSFLAG_JACKET ; Flag that this is a jacket procedure in 
6D Q0000000'GF 9E 0B69 MOVAB G*MTHSSJACKET_HND, (FP) 

0870 ; set handler address to jacket 
are ; handler 
0B70 86795 3; case of an error in special JSB routine 

50 08 BC 70FD 0870 796 MOVH ax(AP), RO ; RO/R3 = arg 

30 0875 797 BSBW MTHSHATAND_R8 ; Call special HATAND routine 
04 BC 50 7DFD MOVO RO, @Arctangent (AP) ; Store result in first argument | 
04 RET ; Return to caller 


- 


SS —_ —y 
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443 o¢ am -SBTTL MTHSHATAND2 = Standard H Floating Point Arctangent With 2 Arguments | 
Bre 04 ; FUNCTIONAL DESCRIPTION: | 
a: 06 HATAND2 =H floating point arctangent function 
OB7E 08 ; HATAND2(X,Y) is computed as following: 
OB7E 810 ; If ¥ = 0 or X/¥ > 2eei14, HATAND2(K.¥) = 90» (sign x) 
OB7E 811 ; If and X/Y =< 2*#1{4, HAT AND G(X, ¥) = HATAND(X/Y) | 
pare o1¢ : If Y < 0 and X/Y¥ =< 2%#114, HATAND2(X,Y) = 180 * (sign X) + HATAND(X/Y) | 
OByE BI ; | 
OB7E 815 ; CALLING SEQUENCE: 
OB7E 816 : | 
OB7E 817 : Arctangent2.wh.v = MTHSHATAND2(x.rh.r, y.rh.r) 
OB7E 818 : | 
OB7E B33 3 ; -or- 
0B7E 20 ; 
pare es! ; CALL MTHSHATAND2(Arctangent2.wh.r, x.rh.r, y.rh.r) 
OB7E 8 : ; Because an H floating result cannot be expressed in 64 bits, it is 
QOB7E 824 ; returned as the first argument, with the input parameter displaced 
OB7E 825 : to the second argument, in accordance with the Procedure Calling 
0B7E Bg8 : 3 Standard. 
OB7E 827 ; 
OB7E 828 ; INPUT PARAMETERS: 
0B7E 4 
00000008 OB7E 830 x = ¢ * LONG ; X is the first argument 
0000000C pare o3 y = 5 * LONG ; Y is the second argument 
pare 635 ; SIDE EFFECTS: See description of MTHSHATAND 
OB7E 835 : OUTPUT PARAMETERS: 
OB7E 836; 
00000004 OB7E 837 Arctangent2 = 1 * LONG ; Arctangent2 is the result 
OB7E 838 ; 
OB7E 839 ;-- 
OB7E 840 
OB7E 841 
41FC OB7E Hk -ENTRY MTHSHATAND2, ACMASK ; STANDdard colt: “ir fererense qeter 
0880 84 ; Disable DV (and FU), enable IV ; | 
pas B44 MTHSFLAG_JACKET ; Flag that this is a jacket procedure in | 
6D 00000000'GF 9E 0880 MOVAB G*MTHSSJACKET_HND, (FP) 
0887 ; set handler address to jacket 
pas? : handler 
pee 845 ; Case_of an error in special JSB routine | 
D 0887 se8 MOVH ax(AP), RO ; RO/R3 = argl | 
D beat Se MOVH @y (AP), R4 ; R4&/R7 = arg2 
iad bcs : Test if Y = 0 or X/¥ > 2e#114 | 
13 0891 «851 * BEQL DEG ; Branch to INF_DEG if Y = 0 
AB 93 3¢ BICW #"x80 0, RO, R8 ; RB = exponent (X) 
AB 99 5 BICW #*x8000, Ra -(SP) ; Stack = exponent(Y) 


ws 
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58 A2 OB9F 54 SUBW2 (SP)+ ; RB = exponent(X) = _exponent(yY) 
0073 8F af Bi OBA 55 CMPW RB rise 3 Compare R with 115 , 
2c 14 BA? 36 BGTR ne _DEG : If K/Y > 2"*114, branch to INF_DEG | 
OBA9 58: Test if Y>Oory<0 
QBA9 3 ; 
54 BS OQOBA9 860 TSTW R4 ; Test the sign of Y 
34 14 QBAB 861 BGTR A2PLUSD 3; Branch to A PLUSD if Y>0 
0 85 QBAD B66 TSTW ; Test the s sien 8 
OE 18 peat yt BGEQ A1PLUSD ; Branch to A PLUSD- if X >= 0 
+} 86> : Y <0 and X < 0 and X/Y =< 2114 
4A 10 0881 867 BSBB “eto _R8D ; RO/R3 = HATAND(X/Y) 
50 FOC8 CF 62FD 0BB3 868 SUBH2 ; RO/R3 = -180 + HATAND(X/Y) 
04 BC 50 7DFD OBB9 869 MOVO RO, aArctangent2(AP) ; Store result in first argument 
04 OBBE 870 RET ; Return to caller 
OBBF 871 ; 
OBBF are 3; Y <0 and X > O and X/Y =< 28114 
OBBF 873 ; | 
OBBF 874 AIPLUSD: 
3C 10 OBBF 875 BSBB or _R8D ; RO/R3 = el 
50 FDBA CF 60FD OBC1 876 ADDH2 ; RO/R3 = 180 + HATAND(X/Y) 
04 BC 50 7DFD OBC7 877 MOVO RO, aArctangent2(AP) ; Store result in first argument 
04 OBCC 878 RET ; Return 
OBCD 879 ; 
OBCD 880 ; Y > O and X/Y =< 2*#114 
OBCD 881 ; 
OBCD 882 A2PLUSD: | 
2E 1 OBCD 883 BSBB MTHSHATAND_R8D ; RO/R3 = HATAND (X/Y) 
04 BC 50 7DFD OBCF 884 MOVO RO, @Arctangent2(AP) ; Store result in first argument | 
0BD4 885 RET 3; Return 
0OBD5 886 ; 
0BD5 887 ; Y = 0 or X/Y > 2e*114 
OBDS 888 ; 
OBD5 889 INF_DEG: 
50 B85 OBDS 890 TSTW RO ; Test the sign of X 
OA 14 QBD7 891 BGTR 1$ : Branch if X > 0 
10 13 OBD9 4 BEQL 2$ : Branch if X = 0 
04 BC FD90 CF 7DFD OBDB 89 MOVO  H_M90, @Arctangent2(AP) Miao | 
OBE2 894 ; Store result in first argument 
04 OBES S32 RET ; Return 
04 BC FD7B CF 7DFD OBES 897 18: MOVO  H_90, @Arctangent2(AP) | 
OBEA 98 ; Store result in first argument 
04 OBEA 899 RET ; Return 
OBEB 900 | 
EB 901 ;+ : 
4 Ons ; Here if both X = 0 and Y = 0. Signal INVALID ARG TO MATH LIBRARY 
EB 904 | 
50 01 OF 79 QBEB 905 28: ASHQ = #15, #1, RO ; RO/R3 = reserved operand, co180ed 
EF 2n8 3 to CHF$_MCH_SAVRO/R1 so handlers | 
BEF 90 ; can change Tf they want to continue. | 
52 7C OBEF 908 CLRQ = R2 : 
7E OO°'8F 9A OQOBF1 909 MOVZBL page + INVARGMAT, -(SP) : Code for invalid argument to math Library | 
00000000'GF 01 FB OBFS 910 CALLS , G*ATHSSSIGNAL : Signal SEVERE error | 


3 
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04 OBFC 911 RET ; Return if a handler says SS$_CONTINUE 
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ATHSHATAND_ = Special HATAND routine 6-SEP-1984 ” 


r 4 | 
MTHRTL.SRCJMTHHATAN. MAR; 1 (12) | 2-00 


Logic for positive medium sized arguments. Get pointer into HATAND_TABLE. 


are $17 «SBTTL MTHSHATAND_R8 = Special HATAND routine 
A B12 ; Special HATAND - used by the sTANDdard routine, and directly. 
BFD 319 ; CALLING SEQUENCE: 
OBFD 318 5 save anything needed in RO:R8 
QOBFD 919 ; MOVH osee © : Input in RO/R3 
OBFD 920 ; JSB MTHSHATAND_R8 
ad 4 ; return with result in RO/R3 | 
OBFD 925 ; Note: This routine is written to avoid causing any wateger overflows, floating | 
yore 4 : 3 ever Tiaet, or floating underflows or divide by 0 conditions, whether enabled or 
3 not. 
OBFD 9 § 3 
OBFD 927 ; REGISTERS ysee 
FD 8928 ; RO/RS = floating argument then result 
D 929; R4/R7 = scratch 
D 930; R8 - Pointer into HATAND_TABLE 
D 931; RO:R5 = POLYH: 
D 336 3 stack - Y during POLYH 
D 93 
D 934 
D 935 
D 936 MTHSHATAND R&D: ; For local use only! 
50 54 66FD > > Bae IVH2 R4, RO 
1 938 MTHSHATAND_R3:: : Special HATAND routine 
50 73FD 939 TSTH RO ; (SP) = X = argument 
3 18 940 BGEQ § POS_ARGD ; 
0081 +331 o03 BRW NEG_ARGD ; Branch to negative argument logic 
308 ; Argument is positive 
944 ; 
945 POS_ARGD: 
Ss fm 8 946 ROTL #3, RO, RB ; RB = Biased exponent bits and first 
58  FFF80000 8F CA 947 BICL #*XFFF80000, RB 3 three fraction bits 
58 OOO1FFEC 8F C2 948 SUBL #°XTFFEC, RB ; Argument is less than 3/32, 
-— 949 BLSS SMA 3 branch to small argument logic 
58 37 «ODT 950 CMPL #°X37, RB ; Argument is greater than of equal to 
4A 15 23) BLEQ LARGE _ARGD 3 11, branch to large argument logic 
i 
57 Q0000000'GF oD 955 MOVAL G“MTHSSAB_ATAN_V, R7 ; R? = Address of RTL vector entr 
57 00000000' GF is 936 ADDL G*MTHSSAB_-ATAN-V, R7  ; R7 = Address of MTH$$AB_ATAN table 
5 6748 95 MOVB (R7)CR8), ; R8 = Index into HATAND_TABLE. 
58 FBA6 CF48 7EFD 938 MOVAO HATAND_TABLECR8], R8 =: R8 = pointer to XHI | 
44 ; Compute z | 
8 306 MOVO (RB)+, R4 3; R4 = XHI 
7E 4 96 SUBH3 = R4, RO, -(SP) > (SP) = X = XHI | 
4 964 MULH2 4, RO ; RO = X*XHI | 
965 ADDH #1, RO ; RO = 1 + X*XHI 
96 DIVH RO, (SP) : (SP) = 2 = (X = XHID/(1 + X*XHI) | 


+ © + + + 6 + 6 + + + > + + > 6 9 + + + +o + ++ + > oo eee TTT TTT Tp ty 
PAOAOOAIOOAOOOAOOOAOOAAOOOOAOOOOOMOAOOCOCBWWWIOIVIOwOaw 


DAS & FAW MWMAIAIWIVRIPPINPINIDEY | — QQOOOOCOOOO Fh “" 
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967 ; 
308 ; Evaluate Z*P(7**2) 
969 ; 
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50 6€ FD C 70 MULH3 (SP), (SP), RO ; RO = Z#*2 
FB24 CF OE $f Serb i er POLYH Ro. MHATANOLEN1=1, led) pe sielsts 
3 = ee 
50 § 64FD C 5 358 MULH (SP)+, RO 3 RO = HATAN(Z) = Z*P(Z*22) 
50 60FD C6 974 ADDH (R8)+, RO ; RO = HATAN_XHI_LO + HATAN(Z) 
50 68 60FD C64 975 ADDH (R8), RO : RO = HATANTX) = HATAN_XHI_HI + 
68 378 : (HATAN_XHI_LO + RAATAR(Z)) 
05 0Cé6 97 RSB 3; Return 
0c69 978 
584 979 
009D 31 O0C69 980 SMALLD: BRW SMALL_ARGD ; Dummy Label used to avoid adding 
Beer 4 1 : on gates ete ag en, in the 
medium argument Logic 
OCee aes an iti t Logi | 
3; Large positive argument logic. 
OC6C 985 ; . . 
OC6C 986 
0céC 987 LARGE_ARGD: 
7E F38E CF 50 67FD Begs 988 DIVH3 RO, H_M1.0, -(SP) ; (SP) = -W = -1/x 
50 6E «46 6E 65FD «0C73)~=—(989 MULH3 (SP), ~(SP). R ; RO = Wee2 
FBO1 CF OF 50 75FD 0C7B8 990 POLYH RO, #HATANDLEN1-1, HATANDTAB1 
OC7F 991 3 RO = P(Wee2) 
50 BE 64FD OC7F 93 MULH2 (SP)+, RO 3 RO = -HATAND(Z) = -Z*P(Wee2) 
50 FCD8 CF 60FD pees 99 ADDH2 4H_90, RO : RO = HATAND(X) = 90 = HATAND(Z) 
05 OC89 994 RSB 3; Return 
OCBA 995 
pega 996 ; 
C8A 997 ; Logic for negative arguments 
OCBA 998; 
OC8A 999 
OC8A 1000 NEG_ARGD: | 
58 50 03 9C QC8A 1001 ROTL #3, RO, R8 ; R8 = Biased exponent bits and first 
58 FFF80000 BF CA OCcBE 100¢ BICL  #XFFF80000, RB : three fraction bits 
58 OOOSFFEC 8F C2 oc 100 SUBL  &#*XSFFEC, RB ; Argument is less than 3/32, 
68 19 OC9C 1004 BLSS SMALL _ARGD : branch to small argument Logic 
58 37 D1 OC9E 1005 CMPL a*x R8 ; Argument is greater than of equal to 
48 15 pCa 1906 BLEQ N_LARGE_ARGD 3 11, branch to Large argument logic 
OCS 1008 ; Logic for negative medium sized arguments. Get index into HATAND_TABLE. 
57 Q0000000'GF ODE OCA 1010 MOVAL G*MTHSSAB_ATAN_V, R7 ; R7 = Address of RTL vector entr 
57 00000000'GF CO OQOCAA 1011 ADDL G“MTHSSAB_ATAN_V, R7 ; R7 = Address of MTH$S$AB_ATAN table 
58 6748 90 OCB1 1 \ MOVB (R7)CR8), ; R8 = Index into HATAND_TABLE. 

58 F825 CF48 7EFD sit ! 1? MOVAO HATAND_TABLECR8), R8 : RB = pointer to XHI | 
CBC 1015 ; Compute Z 
cBC 1318 : 

54 88 7DFD OCBC 101 MOVO.  (RB)+, R4 > RG = XHI 
7E 50 54 61FD cc 1318 ADDH3 = R4, r6, - (SP) ; RO = X + XHI = X = (=XHI) 
50 28 64FD OCC 191 MULH » RO ; RO = X*XHI 
50 08 0 63FD OCC9 1020 SUBH RO, #1, RO 3 RG = 1 = X*XH] = 1 + X*(=xXHI) 
6E 50 66FD CC 1084 DIVH2 RO, (SP) : (SP) = 2 | 
. | 
$8 i 3; Evaluate Z*P(Z**2) 
50 6€ oF a4 CD 1025 MULH3 (SP), (SP), RO ; RO = *2 
FAA2 CF OE SFD OCD? 1026 POLYH RO, #HATANDLEN1-1, HATANDTABI 
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CDE 1027 3: RO = P(Z*#2) 
8 § 64FD CD 1 8 MULH (SP)+, RO 3; RO = HATAND(Z) = Z*P(Z*22) 
eged cE 1 SUBH {ie RO $ = HATAND_XHI_LO + HATAND(Z) 
0 68 62FD CE 1030 SUBH (RB), RO ; RO = HATANDTX) = HATAND_XHI_HI + 
CEA 1031 ; (HATAND_XHI_LO + HATAND(Z)) 
05 ree ! ¢ RSB 3; Return 
CEB 1034 : Logic for Large negative arguments 
CEB 1035 ; . . . . 
CEB 1 § 
CEB 1 N_LARGE_ARGD: 
7E F30F CF 50 67FD OCEB 1 4 DIVHs RO, H_M1.0, =(SP) : (SP) = Wie 1/1X! 
50 6E 6€ 65FD OCFe 1 MULH3 (SP), ~ (SP). : RO = Wee2 
FA82 CF OE 50 75FD CF 1040 POLYH RO, #HATANDLEN1-1, HATANDTAB1 
OCFE 1041 3 RO = P(Wee2) 
50 BE 64FD OCF 1306 MULH (SP)+, RO ; RO = HATAND(W) = WeP(We22) 
50 FCS59 CF 62FD ODO02 104 SUBH H_90, RO ; RO = HATAND(X) = HATAND(W) = 90 
05 4 ee RSB ; Return 
DO9 1045 
0D09 1046 ; 
OD09 1047 ; Small argument logic. 
QOD09 1048 ; 
0D09 1049 
0D09 1050 SMALL_ARGD: 
7E 50 70FD OD09 1051 MOVH RO, -(SP) ; (SP) = argument = X 
2c 13 OD0D 193¢ BEQL 3$ : 
50 8000 8F AA ODOF 105 BICW #*x8000, RO 3 RO = ix! 
50 3FC9 8F B1 OD14 1054 CMPW #*xX3FC9, RO 3; Compare 2*-56 to iX! 
09 19 OQD19 1055 BLSS 1$ ; Needs polynomial evaluation 
50 6E FC30 CF 65FD OD1B 1056 MULH3 H_PI_OV_180_M_64, (SP), R 
Boss 1057 3; RO = X*(pi/180 = 64) 
OF 11 $3 1058 BRB 2$ 3; Argument = 0 = return value 
50 50 64FD OD24 1059 1$ MULH2 RO, RO : = Xee2 
FB41 CF OE 50 75FD O0D28 1060 POLYH RO, #HATANDLEN2-1, HATANDTAB2 
OD2F 1061 : RO = Q(X##2) 
50 6E 64FD OD2F 1388 MULH2 (SP), RO 3 RO = X#Q(X*#*2) 
6 06 AO 0D33 1063 2$ ADDW #*X6, (SP) : (SP) = Xe2e26 
50 8E 60FD 0036 1064 ADDH2 (SP)+, 3 RO = HATAND(X) = X#2e*6 + X#Q(X##2) 
05 OD3A 1065 RSB ; Return 
$038 1066 
8E 73FD OD3B 1067 3$ TSTH (SP)+ ; Remove argument from stack 
05 QOD3E 1068 RSB 3; Return 
OD3F 1069 
OD3F 1070 - END 
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PSECT name Allocation PSECT No. Attributes 
BS . 00000000 0.) O00 ¢ QO.) NOPIC USR CON ABS LCL une NOEXE NORD NOWRT NOVEC BYTE 
MTHSCODE OOOOODSF ¢ 3391.) O01 ¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
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Initialization 29 00:00:00.06 88:80:00 8 
zonsane processing 115 80:00:05: 89 0:00:04.6 

s 1 134 00:00:02.89 00:00:09.91 
Cosel table sort 0 ae SY 00:00:00.04 
Pass 2 215 00:00:02.36 00:00:07.11 
Symbol table output 7 00:00:00.06 00:00:00.18 
Psect synopsis output 3 00:00:00.02 00:00:00.06 
Cross-reference output 0 80:00:06: 45 00:00:00.00 
Assembler run totals 505 00:00:06.10 00:00:22.79 


The working set Limit was 1200 pages. 
19511 bytes (39 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space at tecates to hold 56 non-local and 8 Local symbols. 
1130 source Lines were read in Pass 1, producing 27 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 
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Macro Library name Macros defined 
“$255SDUAZB:(SYSLIBISTARLET.MLB;2 0 (iti‘sé~;~*~*™” 0 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHHATAN/OBJ=OBJ$:MTHHATAN MSRCS:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MSRC 
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